Perturbed disks get shocked. Binary black hole merger effects on accretion disks 
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The merger process of a binary black hole system can have a strong impact on a circumbinary disk. 
In the present work we study the effect of both central mass reduction (due to the energy loss through 
gravitational waves) and a possible black hole recoil (due to asymmetric emission of gravitational 
radiation). For the mass reduction case and recoil directed along the disk's angular momentum, 
oscillations are induced in the disk which then modulate the internal energy and bremsstrahlung 
luminosities. On the other hand, when the recoil direction has a component orthogonal to the disk's 
angular momentum, the disk's dynamics are strongly impacted, giving rise to relativistic shocks. 
The shock heating leaves its signature in our proxies for radiation, the total internal energy and 
bremsstrahlung luminosity. Interestingly, for cases where the kick velocity is below the smallest 
orbital velocity in the disk (a likely scenario in real AGN), we observe a common, characteristic 
pattern in the internal energy of the disk. Variations in kick velocity simply provide a phase offset 
in the characteristic pattern implying that observations of such a signature could yield a measure 
of the kick velocity through electromagnetic signals alone. 

PACS numbers: 



I. INTRODUCTION: 

The study of a number of astrophysical systems will 
soon add gravitational wave astronomy as a new tool to 
complement observations in the electromagnetic band. 
Since most systems capable of producing detectable grav- 
itational waves will also radiate strongly in the electro- 
magnetic band (see, e.g., P, combining information 
from both spectra will allow for a richer description of 
these systems. Furthermore, the complementary nature 
of observation in both bands will help the detection en- 
terprise as a signal in one band will help follow up studies 
in the other (see, for instance, [1,3.) 

Among interesting possible sources of strong signals 
in both spectra, the collision of a binary black hole sys- 
tem within a circumbinary disk presents the possibility 
of a detection of gravitational waves (as the black holes 
merge), which will be followed by electromagnetic sig- 
nals emitted by the disk as its dynamics are affected in 
the process [Bj. This scenario is common in nature, since 
massive black holes exist in the core of most galaxies and 
galaxies undergo mergers throughout their evolutionary 
path. As galaxies merge, they produce a binary black 
hole in the newly formed galaxy which eventually collide 
as their orbit shrinks through several mechanisms. As 
discussed in a circumbinary disk is formed as the 
binary hollows out the surrounding gas, and the disk 
becomes mostly disconnected from the binary's dynam- 
ics 0,01 • Afterwards, while the disk remains essentially 
frozen, the black holes' orbits continue to shrink until 
they merge. 

The merger process gives rise, in particular, to two 



relevant effects that will perturb the disk (see, e.g., [1, 
@). One is related to the final mass of the black hole, 
which is less than the initial total mass as the system 
radiates energy via gravitational waves^. The other one 
is a consequence of the radiation of linear momentum, 
which if asymmetric (as in the case of an unequal mass 
binary, or asymmetric individual angular momenta of the 
black holes), induces a nontrivial recoil on the nascent 
black hole. This recoil effect has been predicted before 
through perturbative analysis of Einstein equations [lol . 
ITl| . and recent numerical simulations implementing the 
equations in full show even higher recoil velocities are 
possible [H, [H, M, m, m, 113, M, [H, HO, mi . The largest 
recoil velocities found correspond to mass ratios close to 
1 and spins lying anti-aligned on the orbital plane. In the 
case of quasicircular orbits, recoil velocities up to about 
4000 km/s have been calculated [ij. However, most of 
the black hole collisions occurring in nature are expected 
to produce kicks of about 500 km/s or less, since larger 
kicks would occur only in the case of nearly equal masses 

As a result of both effects mentioned above, the fluid 
dynamics in the disk is modified and shocks may be in- 
duced. The shocks' energy can then heat the gas, which 
can produce electromagnetic flares. These flares are ex- 
pected to occur later (a few months to years) , and to last 
considerably longer (thousands to hundreds of thousands 



^ Possible observable consequences of this effect were first dis- 
cussed in 01 . 
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of years) [1^ [U, [2^ , suggesting tantalizing prospects for 
LISA observations aiding and complementing the electro- 
magnetic observational prospects of these systems. For 
the recoiling black hole case, prior studies, which employ 
simulations of coUisionless particles in Keplerian orbits 
forming a flat (zero hei ght] disk, predict emissions rang- 
ing from UV to x-rays [2^, [13] or in the infrared [1^ if 
this radiation is assumed to be absorbed before leaving 
the disk and re-emitted. Since these studies employ a 
particle description of the fluid, they can not fully cap- 
ture the development (and hence influence) of shocks, 
which must be estimated by detecting collisions between 
particles. A recent work [2g adopted a field description 
for the fluid and studied the impact of a mass reduction 
in a pseudo-Newtonian potential to account for an in- 
nermost stable circular orbit (ISCO) at r = 6M (which 
corresponds to the ISCO of a nonspinning black hole, 
while this is a rather uncommon output [27l . [28l . [29l . [soj 
in the merger of two black holes, the spin value will play 
a relevant role mainly if accretion develops). Based on 
computations of bremsstrahlung luminosity, that work 
predicts a decrease in luminosity as the fluid orbits ad- 
just to the reduced gravitational potential. 

In this work we study the effects on the disk by also 
considering a perfect fluid but in our case we do so em- 
ploying the fully relativistic hydrodynamic equations in a 
background space-time. Thus, we are able to examine ef- 
fects of spin, mass reduction and accretion, and comment 
on the relevance of different processes. In particular, our 
studies indicate that a significant distortion of the disk 
develops as time progresses when the kick has a compo- 
nent perpendicular to the disk's axis and that qualita- 
tively similar features are present in all these cases. 

In Section II, we briefly review our formulation of the 
problem and numerical approach. Section III describes 
our initial configuration. We discuss the observed dy- 
namics in the disk after the merger has taken place in 
section IV, taking into consideration the effect of mass 
reduction and different recoil velocities. Section V con- 
cludes and offers some further considerations. 

Unless otherwise specified, we use geometrized units, 
where G ~ c — 1, and sum over repeated indexes. Greek- 
letter indices range from to 3. 



II. OVERVIEW OF THE NUMERICAL 
APPROACH: 

We implement the general relativistic MHD equations 
using a hig h resolution shock capturing module described 
in [3l|, [33 . We introduce however a slight modification 
of the hydrodynamic equations inside the horizon to im- 
prove the fluid's behavior close to the excision region. 



Given an equation of the form^ 

ij + F{uy = s, (1) 

we modify it in order to include a damping term: 

U + F{Uy = S-f{r){AxnU-Uo), (2) 

where the function f(r) decreases smoothly with r, from 
100 at the excision region to zero at the event horizon 
(EH), and is zero for r > rEH, so that the exterior of the 
BH is causally disconnected from the effect of this extra 
term. Uq is set to zero or to the value of the atmosphere 
if the corresponding field has one. The coefficient (Ax)^ 
ensures that the damping term converges to zero and will 
not modify the convergence rate as long as one chooses 
p to be greater than or equal to the order of convergence 
of the code. In this work we adopt p ~ 4. 

These equations are implemented within the had 
computational infrastructure which provides distributed 
Berger-Oliger style adaptive mesh refinement (AMR) [H, 
[S^ l with full subcycling in time, together with a novel 
treatment of artificial boundaries [11] . Because of the dy- 
namics involved in this work, it is only necessary to use a 
fixed refinement hierarchy, covering with finer grids the 
(central) region containing the disk and black hole, and 
increasingly coarser grids in the outer regions in order to 
locate the boundaries far away at a low computational 
cost. 



HI. OVERVIEW OF THE PHYSICAL SETUP: 

To explore the effects of the black hole merger in the 
dynamics of the accretion disk, we concentrate, in partic- 
ular, in the post-merger stage -when the main burst of 
gravitational radiation has passed through the disk and 
this has settled down to a quasistationary state-. To 
simulate a BH formed through the merger process and 
account for the main effects of mass loss or recoil, we ei- 
ther consider a reduction in the mass of the black hole by 
5% or apply a boost to the BH in a given direction. In 
the latter case, it is easier to adopt the BH's rest frame 
and apply the boost to the fluid variables (in the oppo- 
site direction) describing the disk, which is represented 
by a stationary toroidal solution of the fluid equations 
in a Kerr background. Thus, starting with a stationary 
torus on a Kerr background, we perform a Lorentz boost 
with velocity — Wkick on the disk. We employ this boost 
to transform the fluid's 4-velocity and magnetic field 
4-vector b'^ when considering the recoil case. 

The toroidal solutions are constructed following an ap- 
proach similar to that in [sl], adapted to the ingoing 



^ For clarity we use here a simple expression to represent the fluid 

equations. See |31| . |32|| for the full equations. 
^ Studies of possible premerger effects are presented in. e.g., [ssl . 

m. 
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Kerr-Schild coordinates adopted in our studies, and with 
a different choice of specific angular momentum for the 
fluid for easier comparison with previous work in the ab- 
sence of magnetic fields . While in the current work we 
do not simulate scenarios that include a magnetic field, 
we discuss the construction of initial data that allows 
for doing so for future reference. In our case, we adopt 
the more standard I = —u^ iut)~^ = const, (see below) 
to allow for an easier comparison with previous work in 
the absence of magnetic fields. In particular, we verify 
that identical solutions to those of [331 obtained if the 
magnetic field is set to zero. In what follows we review 
the main steps in this construction. 

The stress-energy tensor for ideal MHD can be written 

as 



T, 



(3) 



where p, P, and are the fluid's density, pressure and 
4-velocity, respectively, h^^ is the magnetic field 4-vector, 
and h is the specific enthalpy, defined as 



p/i = p(l + e) + P, 



(4) 



where e is the specific internal energy density. 

For the construction of initial data, we work with cylin- 
drical coordinates (t, r, <j), z) and make the assumption 
that the space-time is stationary and axially symmet- 
ric. We adopt coordinates adapted to these symmetries, 
so that only the t- and ^-components of u^^ and are 
nonzero. 

The fluid equations arc obtained from 



(5) 



together with the continuity equation Vp(pw^) = 
(which is trivially satisfled under our assumptions). Af- 
ter some manipulation, equation ([5]) can be reduced to 
the integral equation 



hp 



2phD 



where D = \gttg<p4>—gtcji\- This equation can be integrated 
after imposing further conditions that flx relationships 
between the fluid variables as discussed below. 

First, we fix a relationship between the velocity com- 
ponents. This can be accomplished by requiring that the 
specific angular momentum / satisfles 



I 



U4, 

Ut 



const. 



(7) 



Second, we assume an isentropic fluid, imposing dh ~ 
p^^dP, which allows us to integrate one of the terms out. 
An equation of state that satisfies this condition is that 
of a polytrope 



(8) 



In this case, the specific internal energy density can be 
calculated as 



r - 1' 



(9) 



We adopt this condition only to obtain the stationary 
solutions for initial data. The fluid's entropy will change 
after the kick and so we adopt, during the evolution, a 
F-law equation of state 



P= (F-l)pe 



(10) 



with F = 5/3 considering the gas as being monoatomic. 

Finally, we impose a convenient expression for 6^ in 
terms of other variables to integrate the last term 



I? D = C'iphD) 



(11) 



where C and q > 1 are arbitrary constants. 

After integrating eqn. (6), we use ([8]) and ([9]) to ehm- 
inate p and e and obtain an algebraic equation for P, of 
the form 



F{P,g^,,l,C,q) = Fo 



(12) 



where Fq is a constant of integration. This equation can 
be solved analytically in the absence of magnetic field 
{b^ = 0), otherwise a straightforward numerical integra- 
tion can be set up to obtain the solution. The boundary 
of the torus is determined by setting P = 0, obtaining an 
expression of the form 



(13) 



d(b D) = const., (6) condition u^'u 



which, through the dependence of g^^, on the coordinates, 
is an implicit surface equation. Notice that it is indepen- 
dent of both C and q so that the location of the disk's 
boundary is independent of the magnetic field. The so- 
lutions obtained may be toroidal as well as spheroidal, 
depending on the values of I and Pg- 

Once P is known, one can use once again equations 
(|H1) and © to recover p and e. The velocity is ob- 
tained from equation ^ together with the normalization 



-1. Finally, the magnetic field 6^ is 



determined by eq uation pT|) together with the relation 
Wu^ = (see [311 )■ The magnetic filed is always zero at 
the surface of the disk (from equation [TT]), and one can 
control how rapidly it decays to zero with the parameter 
q, and its maximum magnitude with C. 

As mentioned, when considering the recoiling case, 
the initial data for a disk is given a Lorentz boost 
with respect to the stationary system of the background 
black hole. A representative example of the toroidal 
configurations is shown in Figure [TJ Such configu- 
rations are then evolved on a computational domain 
given by [— 150M, 150 M] in the x- and y-direction, and 
[—100 M, 100 M] in the z-direction (since the disk lies on 
the x — y plane), with an FMR configuration having 3 lev- 
els of refinement. The code used in our studies has been 
previously tested and employed in a variety of stringent 
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scenarios, e.g.jjH, [13, |4l|. For our specific application 
we have verified that in the absence of a kick or mass re- 
duction the disk remains stationary during the evolution 
as expected. Additionally, we have verified convergence 
by comparing results obtained with three different reso- 
lutions in the case of a kick velocity of 3000 km/s per- 
pendicular to the axis of rotation. The convergence rate 
measured at different locations varies between first and 
third order depending on the presence of shocks. 

Certainly the parameter space is too vast to allow an 
exhaustive computational study. Therefore, wc mainly 
concentrate here on varying the most relevant parame- 
ters, i.e., the kick magnitude and direction and study a 
few other cases varying the spin parameter to verify our 
results are qualitatively the same. Notice that variations 
with respect to the spin parameter a should not lead 
to significant qualitative differences unless accretion de- 
velops, as the disk's inner edge is located sufficiently far 
away for its influence to be of higher order. This intuitive 
observation is confirmed by our simulations. 



bound to the black hole in all cases considered^. In 
fact, the binding energy per unit mass at the surface of 
the torus is 0.0121c^, which implies a escape velocity of 
0.155c. 

Notice that the location of the disk's inner radius 
can vary significantly depending on diverse physical 
parameters (e.g. kinematic viscosity of the gas, accretion 
rate, binary mass ratio, etc.) We adopt a small value 
but within the allowed ones to reduce the computational 
cost of the long simulations required and concentrate on 
extracting physically robust conclusions, which can be 
intuitively extended to general cases. 

Throughout the rest of this paper, unless otherwise 
specified, all kick orientations mentioned refer to the kick 
(or Lorentz boost) applied to the disk, which would cor- 
respond to the black hole being kicked in the opposite 
direction. 



IV. RESULTS: 
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FIG. 1: Representative example of the toroidal initial con- 
figurations, showing the density at the equatorial plane (left 
panel), and at a meridional plane (right panel). The dashed 
line indicates the location of the event horizon. 



The toroidal solutions employed in this work corre- 
spond to specific angular momentum l/M ~ 6, spin 
parameter a/M = 0.5 (except when analyzing the 
solution's dependence on the spin where we also consider 
a/M = 0.9). Also, wc fix the magnetic field parameter 
C = (so that 6 = 0), and choose Fq so that the 
inner edge of the disk is located at rjn = 20M. With 
this choice of parameters, the outer edge is located at 
Tout = 60-Af and the maximum pressure in the disk 
lies at 7'in = 33M. The orbital velocity of the fluid is 
then 0.28,0.17 and 0.10c at rin,rm and rout respectively. 
Thus, the orbital period at fm is Pm = 1220M. The 
sound speed has a maximum value < 0.05c close to r,n, 
and drops abruptly to zero at the boundary of the torus. 
All fluid elements in the torus have an orbital speed 
much greater than the highest kick velocity adopted 
in this work, i.e., 0.01c=3000 km/s and so will remain 



A. Diagnostic quantities: 

We monitor the fluid's behavior by examining the de- 
pendence of the primitive values as different physical pa- 
rameters are varied. Ultimately, our goal is to understand 
possible electromagnetic signals emitted by the system 
as the disk' dynamics is affected. At present, our simula- 
tions do not incorporate radiation transport; thus, a di- 
rect computation of these signals is not possible. There- 
fore, we concentrate on related quantities, which when 
combined with a suitable model, can be tied to possible 
emissions. In particular, we compute (an approximation 
to the) temperature (T), the total internal energy (U) 
and bremsstrahhmg luminosity {Lb) as 

T ^ P/p, (14) 
[/ cx y pedV, (15) 

Lb oc J p^T^/^dV . (16) 

Notice that unless the disk is optically thin, the 
bremsstrahlung luminosity need not capture the lumi- 
nosity resulting from shocks and shock heating. While 
the bremsstrahlung luminosity is a good measure of the 
energy exchanged between atoms and the radiation field, 
it does not take into account how this energy can be ra- 
diated out of the disk. In the absence of a more refined 
model, the qualitative features of the true radiative be- 
havior can be estimated simply by a black body assump- 
tion. We thus monitor the internal energy for this pur- 
pose and also the bremsstrahhmg luminosity to obtain a 



* For comparison purposes we have also employed the unrealistic 
value of 9000 km/s 
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measure of the mentioned energy exchange (as well as to 
make contact with results presented in [26|). 



B. Axisymmetric cases: Black hole mass loss and 
kick along disk's angular momentum 

As a first step we consider the effect of BH mass 
loss and that of a kick along the disk's orbital angular 
momentum. The former entails solely decreasing the 
mass of the black hole, while for the latter the mass is 
unchanged but a kick is introduced along the z axis. In 
both cases, the underlying axisymmetry of the problem 
is not broken, which as we shall see later, is a key issue. 

For the mass loss case, we employ a toroidal solution 
corresponding to a black hole of mass Mq for the ini- 
tial data, and set M = O.OSAfg. The dynamics of the 
disk with either a reduced mass or a kick along the z 
axis behave in a rather smooth manner. For the mass 
reduced case, radial oscillations are induced as the differ- 
ent fluid elements follow their corresponding epicycles. 
For the case with a recoil velocity, further oscillations 
are generated by induced motions in the z axis as il- 
lustrated in figure [H Indeed, the recoil motion of the 
black hole introduces a time-dependent vertical compo- 
nent of the black hole's gravitational pull on the disk. 
Using Newtonian mechanics for simplicity and ignoring 
pressure forces, one can show that a particle on a circular 
orbit with velocity Worb, after a vertical kick of magnitude 
I'kick only reaches a height z — '\/2-R(vkick/i^oib) above the 
original plane before turning around. Since Wkick is the 
same for all disk radii, the vertical displacement is mini- 
mal at rjn and maximal at rout- This results in a flexing 
axisymmetric mode, with the outer edge flopping about 
the most. This is supported by Fig. [2] if one defines the 
"midplane" of the disk by joining points at which the 
contours are vertical. This argument ignores pressure, 
but pressure gradients are unlikely to be very important 
away from shocks, and the behavior is qualitatively the 
same. Because all particles on a given radial annulus are 
kicked simultaneously, they remain in phase with each 
other and the flexing mode is naturally excited. Note 
that maximum compression occurs twice per orbital pe- 
riod so this flexing mode is visible in both internal energy 
and bremsstrahlung at a frequency of about twice the or- 
bital (See Figs. [3] and El). 

Most importantly, in either case no significant shocks 
are developed during the time of these simulations (~ 
6Pni )■ The observed smooth behavior translates into 
a rather monotonic behavior in our diagnostic vari- 
ables. Figures [3] and m illustrate the internal energy and 
bremsstrahlung luminosity, respectively. The behavior 
observed in the latter case is qualitatively similar to re- 
sults shown in [2^, i.e., an initial drop followed by a 
recovery in luminosity. Our simulations, which extend 
farther, indicate that this behavior continues quasiperi- 
odically. Notice however that the disk geometry consid- 



ered in [26| is different from ours and the bremsstrahlung 
computed includes the inner portion of the disk while we 
do so for the complete disk. Last, the small drift observed 
in figure [3] is consistent with a linear accumulation of nu- 
merical error. A similar linear drift is observed in sim- 
ulations of an unperturbed disk. This growth however 
is small -within 5%- over the length of the simulations 
considered (~ 6Pni ) and significantly smaller than the 
effects induced by the perturbations due to the recoiling 
black hole. 
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FIG. 2: Density at plane y = in the case of a disk kicked with 
a velocity 3000 km/s in the positive z-direction. The panels 
show snaps from tjM = 500 (top left) to 2000 (bottom right) 
at Af/Af = 500 intervals. 



C. Asymmetric cases: Kick with component 
orthogonal to disk's angular momentum 

Next we concentrate on the oblique recoil case. For 
concreteness we adopt recoil velocity values Vkick = 300, 
1000 and 3000 km/s (we also consider 9000 km/s to verify 
the appearance of the main feature and check the empir- 
ical law presented below). We begin by examining the 
case where the kick direction is on the orbital plane (i.e., 
orthogonal to the angular momentum of the disk). The 
simulations for the different cases proceed along qualita- 
tively similar phases, which are illustrated for the case 
of Wkick=3000 km/s in figure for p at z = 0, and in 
figure [H] for |VP| at 2; = 0. The asymmetry introduced 
by the kick's direction induces an accumulation of gas at 
one side of the disk, while causing a significant decrease 
on the opposite side. As time progresses, shocks develop 
and a complex dynamic arises, at late times ~ 6000Af , an 
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t/M 

FIG. 3: Total internal (normalized) energy of the disk. The 
continuous line corresponds to a BH mass loss of 5% and no 
kick, while the dashed line corresponds to a kick with velocity 
Uk = 3000 km/s along the axis of rotation (and no BH mass 
loss). The vertical scale and range was chosen to coincide 
with those in the other energy plots in this work for easy 
comparison. 




FIG. 4: Bremsstrahlung luminosity (normalized) of the disk. 
The continuous line corresponds to a BH mass loss of 5% 
and no kick, while the dashed line corresponds to a kick with 
velocity Hk = 3000 km/s along the axis of rotation (and no 
BH mass loss). 



accretion phase is clearly noticeable for Wkick >1000 km/s 
(see figure [71). 

To analyze the impact of the disk dynamics and pos- 
sible observable features, we compute the internal en- 
ergy (Fig. ^ and bremsstrahlung luminosity (Fig. [TU)) 
for Wkick=300, 1000, and 3000 km/s. An initial relatively 
small bump is observed, which takes place at a time given 
by half the orbital period of the maximum density re- 
gion, which is consistent with the epicyclic picture. From 
there on, a complex behavior is observed, though no- 
tably, irrespective of the magnitude of the kick, the same 
qualitative features are observed -especially in Fig. [S] 
Generally, we see that both the internal energy and the 
bremsstrahlung luminosity dip and rebound but the in- 




FIG. 5: Density p at plane z — in the case of a disk kicked 
with a velocity 3000 km/s in the positive x-direction, i.e., to 
the right of this page (which corresponds to the black hole 
being kicked to the left). The panels show snaps from t/M = 
500 (top left) to 4000 (bottom right) at At/M = 500 intervals. 



ternal energy ends up higher, while the bremsstrahlung 
luminosity finishes lower. This can be understood as fol- 
lows: the kick energy is dissipated in shocks, increas- 
ing the temperature and the pressure but the subsequent 
expansion reduces the density below the initial values. 
Because the bremsstrahlung emissivity is cx 

p2ri/2, the 

net effect is a reduction in emissivity despite the increase 
in pressure. The relative changes in both internal en- 
ergy and bremsstrahlung luminosity are relatively mod- 
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FIG. 7: Later stages of the simulation shown in figure O in 
which the gas begins to accrete into the black hole. The panels 
show snaps from t/M = 6000 (top left) to 7500 (bottom right) 
at M/M = 500 intervals. Notice that ai t/M = 7000 the 
ISCO is clearly noticeable. 




FIG. 6; To illustrate the formation of shocks we show here 
|VP| at z = 0. Kick velocity of 3000 km/s in the positive 
x-direction. 



est, at a level of 20 — 40% and occur on characteristic 
timescales on the order of lOOOAf = SOOOA'/g s, where A/g 
is the mass of the black hole in 10^ A/0. 

Second, we examine the dependence on kicks at dif- 
ferent angles. Since the main quahtative features of all 
kick cases considered are similar, we concentrate in the 
case ?;kick=3000 km/s as this is the one that displays the 
overall behavior within the shortest computational time. 
We compute the internal energy and bremsstrahlung lu- 
minosity for kicks &i 9 ^ 0,30,60, and 90°, where the 
angle 6 is measured with respect to the axis of rotation. 
Figure [5] shows the density at plane z = in the case of 
a kick at 30° . Figures [TT] and [12] illustrate the (normal- 
ized) internal energy and bremsstrahlung luminosity vs 
time for the different angles considered. Recall that no 
significant shocks form when the kick is along the axis of 
the disk. When the kick has a component along the disk's 
plane however, the qualitative features observed in the in- 
ternal energy are similar for all cases. We note that the 
evolution we observe for a given Vi_ — (3000 km/s) sin ^^ 




-60 ^ -20 20 40 
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FIG. 8: Density p at plane 2: = in the case of a disk kicked 
with a velocity 3000 km/s at 6* = 30°. The panels show 
snaps from t/M = 2500 (top left) to 4000 (bottom right) at 
At/Af = 500 intervals. 



is nicely bracketed by evolutions with pure orbital plane 
kicks above and below Vi_. Thus, dj^ is the most impor- 



8 



1.3- 
1.2- 



„v^= 300km/s 
....V|.= 1000 km/s 
^ -, V, = 3000 km/s 
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FIG. 9: Normalized internal energy. Kicks perpendicular to 
axis of rotation {9 — 90°). 



— vk = 


300Km/s 


— vk = 


lOOOKm/s 


— vk = 


3000Km/s 




FIG. 10: Normalized brenisstrahlung luminosity. Kicks per- 
pendicular to axis of rotation {9 — 90°). 



tant parameter determining the behavior of the kicked 
disk, apart from the small oscillations also present when 
the kick is parallel to the axis of rotation, and the likely 
small differences in the shape of the initial shock. 

Another feature common to all the internal energy (or 
pressure) results (See Figs. |9] and [TT]) is a rapid swing 
from a dip to a bump, followed by an oscillating growth 
at a moderate pace. While the magnitude of the upward 
swing of the internal energy does not depend strongly 
on fj^, the time at which it occurs does. The delay we 
observe decreases as v± is increased. If this delay were 
due to the time taken by a perturbation traveling at v± 
to cross some fixed distance, one would expect a depen- 
dence cx uj^. Instead, we observe a logarithmic decrease. 
Defining the delay as the time after the initial kick at 
which the internal energy swings upward through the ini- 
tial value, we find the following empirical dependence: 



5200 - 912 In 



(17) 



M " ""V 300 km/s 

Note that this formula applied naively "predicts" an in 



finite delay for a kick along the axis of rotation. 

As is well known, constant specific angular momentum 
tori are prone to a nonaxisymmetric corotation instability 
[4^ . I43I I whose nonlinear development has been explored 
numerically in the pseudo-Newtonian approximation (4^ 
and in GRMHD [4^. The final outcome depends on the 
aspect ratio of the torus, the nature and strength of any 
large-scale magnetic fields present, the presence of ac- 
cretion [4^, and this remains to be fully investigated in 
GRMHD context. Therefore, any substantial perpendic- 
ular component of the kick is likely to excite at some 
level the m = 1 nonaxisymmetric mode, which is ex- 
pected to grow at a rate uj « 0.2rini, where rim = 27r/Pni 
is the Keplerian angular frequency at the pressure maxi- 
mum. For the parameters of the torus of our simulations 
Pm = 1220M, and fin, = 0.00515. The behavior de- 
scribed above is suggestive: if one assumes that the ini- 
tial pressure perturbation is 5Pq cx u^. which is reason- 
able for shocks and on dimensional grounds, and one sets 
5P = SPq exp Lut, then the time required for the perturba- 
tion to attain a given fiducial level would follow an equa- 
tion of the form p?|) . with t — t^ci — (2/a;) In (z;_L/z;i.of ), 
where fj-ef and Wrof are some arbitrary reference val- 
ues. Analyzing the results we obtained indicates that 
u! = 0.43il„i, which is on the order of the expected fre- 
quency but significantly higher. Thus, we suggest tenta- 
tively that the swing we see in both the internal energy 
and bremsstrahlung plots in all cases where there is a 
nonzero v± is a common transient response to the kick 
that may be observable in principle, and that the subse- 
quent growth may be due to the growth of the instability 
and/or the rise to the expected level of dissipation of the 
input kinetic energy. At late times for the higher kicks 
our simulations display an accretion phase and so this 
possible saturation can not be explored, though a sugges- 
tive behavior consistent with this saturation is displayed 
by the largest kick considered. 

We note that a similar swing in the bremsstrahlung lu- 
minosity was observed in the (axisymmetry preserving) 
simulations by O'Neill et al [26[ using thin disks, which 
are not prone to the Papaloizou and Pringle instability. 
In the near future, to further elucidate the relative im- 
portance of the transient response and the instability, we 
are planning an investigation of the effects of kicks in 
tori with flatter rotation laws f2 oc r^"^ since the afore- 
mentioned instability does not occur if g < a/S as well as 
examining magnetized tori. 



D. Dependence on black hole spin 

Finally, we investigate possible differences between 
cases with different black hole spins by performing a sim- 
ulation with spin a/M = 0.9 in addition to the value 
a/M = 0.5 used in the rest of the simulations. For this 
test we chose the setting with kick velocity of 3000 km/s 
perpendicular to the disk's axis. Notice that although all 
other parameters coincide in these simulations, the sta- 
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FIG. 11: Internal energy. Kicks at varying inclinations 6 
with respect to the axis of rotation. All cases with Wkick = 
3000 km/s. 
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FIG. 12: Bremsstrahlung luminosity for kicks at varying in- 
clinations 6 with respect to the axis of rotation. All cases 
with tikick = 3000 km/s. 



tionary disk solutions used to construct the initial data 
are slightly different since they depend on a. Still, we 
see no significant differences, as is illustrated in figure [Ql 
where we show a comparison between the maxima of den- 
sity, normalized by dividing by its value at < = 0, which 
is slightly different in each case. 



V. CONCLUSION 

In the current work we have studied the possibility 
that binary black hole mergers, within a circumbinary 
disk, give rise to scenarios likely to emit electromagnetic 
radiation. We have studied both the impact of mass loss 
in the system and recoil velocities. While both induce de- 
formations of the disk, it is the case of a recoiling black 



hole, when the recoil's direction has a component along 
the disk plane, that appears as the most promising op- 
tion to generate an observable, electromagnetic, si gna- 




1000 



2000 3000 4000 
t/M 



5000 



FIG. 13: Comparison of the maximum of density for black 
hole spins a/M — 0.5 (dashed line) and 0.9 (continuous line) 
for a kick of 3000 km/s perpendicular to the disk's axis. 



turc. This is so not just because the effect is larger, 
but also the variability induced is significantly more pro- 
nounced than that observed in the case of mass loss or 
kick along the disk's angular momentum. Furthermore, 
we find that the magnitude of the kick is not very im- 
portant as long as it is less than the smallest orbital 
speed of the fluid. While the kick magnitude impacts the 
time at which the strongest variation in internal energy 
or bremsstrahlung appears, the intensity and time scale 
of the variation and behavior afterwards is not. Since 
supermassive binary black hole mergers will generically 
give rise to recoils (simply by having a mass ratio differ- 
ent from unity) which in turn ensures a kick component 
orthogonal to the final black hole spin, effects like those 
observed here indicate a possible a common behavior for 
the majority of scenarios. 

Our studies also indicate that the final black hole spin 
has no strong effect on the main features of the solution. 
However, if an accretion phase takes place, the location 
of the innermost stable circular orbit will naturally play 
a key role. 
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